Impurity induced spin gap asymmetry in nanoscale graphene 
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We propose a unique way to control both bandgap and the magnetic properties of nanoscale 
graphene, which might prove highly beneficial for application in nanoelectronic and spintronic de- 
vices. We have shown that chemical doping by nitrogen along a single zigzag edge breaks the 
sublattice symmetry of graphene. This leads to the opening of a gap and a shift of the molecular 
orbitals localized on the doped edge in such a way that the spin gap asymmetry, which can lead 
to half-metallicity under certain conditions, is obtained. The spin-selective behavior of graphene 
and tunable spin gaps help us to obtain semiconductor diode-like current-voltage characteristics, 
where the current flowing in one direction is preferred over the other. The doping in the middle 
of the graphene layer results in an impurity level between the HOMO and LUMO orbitals of pure 
graphene (again, much like in semiconductor systems) localized on the zigzag edges thus decreas- 
ing the bandgap and adding unpaired electrons, and this can also be used to control graphene 
conductivity. 



INTRODUCTION 



Applications of graphene with its unique physical 
properties^'^i^'^ in nanoelectronics^'^, magnetism and 
spintronics^iiSiiiiii^, hang crucially on its bandgap and 
spin ordering at the zigzag edges. A bandgap can be 
opened in graphene by breaking the certain symmetries. 
For example, interaction of graphene with its substrate, 
such as SiC, leads to the charge exchange between them 
which breaks the sublattice symmetryi^. Moreover, the 
quantum confinement effect also has been found to intro- 
duce a small bandgap in graphene nanoribbonsii, just as 
was predicted earlier theoreticall y-'^^d^d'^d^ . The effect of 
bandgap opening and spin ordering between the zigzag 
edges are found to be directly linked to each other—. 
When the spins align along the zigzag edges and spin 
states localized at opposite edges have the same spin ori- 
entation, then symmetry of graphene is preserved and 
the system is gapless. Otherwise, if the spin-up states 
are localized along one zigzag edge and the spin-down 
along the other, the sublattice symmetry is broken which 
leads to a gap. In the light of a recent breakthrough in 
fabrication of nanoribbons of required size through unzip- 
ping of carbon nanotubes, the nanoribbons and nanoscale 
graphene are the most promising systems for application 
in nanoelectronics^S. 

Manipulation of the spin ordering is important for 
both graphene magnetism and its electronic proper- 
ties. There are several approaches which have been pro- 
posed recently to control the spin ordering along the 
edges2ii2^i2^i2ii2^i2£i22. One of them is the termination 
of the zigzag edges by functional groups^^. This has the 
advantage that one can achieve half metallicity in this 
process. However, there are some serious issues involved 
here. Firstly, many of these functional groups are placed 
out of the graphene plane thus making the whole struc- 
ture non-planar and, most importantly, termination was 
applied to every second edge cell, which makes its tech- 
nological application very difficult. In fact, we found that 
the strong interactions of the graphene lattice with some 



of the functional groups, such as NII2 and NO2, lead to 
buckling of the graphene layer and twisting of the func- 
tional groups, which subsequently may result in the dis- 
appearance of the half-metallicity of graphene. The curl- 
ing of the graphene layer has been seen in earlier studies 
as well^l, where the boundary conditions were found to 
control the graphene planarity, namely the curling occurs 
for stand-alone systems. For nanoscale graphene the fer- 
romagnetic ordering of the spin states along the zigzag 
edges can be also achieved as subsequence of adsorption 
of gas and water molecules on the graphene surface, as 
we have shown in our previous study22,. The adsorption 
leads to pushing of the a- and /?-spin states to the oppo- 
site zigzag edges thereby breaking sublattice symmetry 
and opening a gap. In some cases the spin asymmetry can 
occur. For example, the adsorption of HF gas molecule 
provide the HOMO-LUMO gap of 2.1 eV for a-spin state 
and of 1.2 eV for /3-spin state. However, due to the weak 
interaction between adsorbant and graphene surface the 
phenomena of the spin alignment along the edges takes 
place locally, thus limiting its application. 

The connection of the phenomena of bandgap opening 
and of the spin ordering with the sublattice symmetry 
lead us to conclude that breaking of this symmetry is the 
main direction to achieve the required semiconductor- 
type bandgap in graphene and a tunable spin ordering. 
Here we make a proposal that the symmetry breaking can 
be done by chemical doping along a single zigzag edge. 
This method is far superior to the earlier approaches in- 
volving edge termination by functional groups because 
doping can be done for every unit cell along the zigzag 
edge and thus preserve the planarity of graphene. Dop- 
ing not only breaks the graphene symmetry, but also can 
induce the spin gap asymmetry due to the energetic shift 
of the molecular orbitals localized on the doped edge. 
In a structure with broken symmetry, the HOMOq and 
LUMO/3 orbital states are localized at one edge, while 
HOMO/5 and LUMOq are at the other. Suppose the 
doping shifts the HOMOq and LUM0;3 orbital states 
localized at one edge down, then the HOMOq-LUMOc 
bandgap (Aq) is increased, while A/j, in contrast, will 



be reduced. If a certain type of impurities can cause a 
significant shift of the bands, then the half-metalhcity of 
graphene may occur. This is what we set out to inves- 
tigate here. An important advantage of this approach 
is that we expect insensitivity of spin selective behav- 
ior to the quahty of the edges, when the band shift in- 
duced by the impurities is stronger than the contribu- 
tion from the edge defects. We also investigate the pos- 
sibility of obtaining an impurity level in the middle of 
the graphene bandgap by doping (in analogy to semicon- 
ductors), which has a lot of technological implications 
as well. Our study of nanoscale graphene was based on 
the quantum chemistry methods using the spin-polarized 
density functional theory with the semilocal gradient 
corrected functional (UB3LYP/6-31G) performed in the 
Jaguar 7.5 program^^. 



II. SYMMETRY OF NANOSCALE GRAPHENE 

Bulk graphene has hexagonal symmetry, while the 
highest possible symmetry of nanoscale graphene would 
be the D2h planar symmetry with an inversion cen- 
ter. The D2h symmetry results in structurally identi- 
cal corners exhibiting ferromagnetic ordering of the spin- 
polarized states localized at the corners, as presented in 
Fig. [1] (a). According to the NBO (natural bond or- 
bital) analysis, the localized electrons at the corners are 
unpaired sp electrons belonging to non-bonded orbitals, 
which are located at the bottom of conduction band or 
top of the valence band. For this symmetry, both a— and 
/3-spin states of the HOMO and LUMO orbitals are local- 
ized on the zigzag edges but their spin density is equally 
distributed between two edges. For nanoscale graphene 
of D2h symmetry the HOMO-LUMO gap appears due to 
confinement and edge effectsi^. The degeneracy of the 
a- and /3-spin states belonging to the HOMO and LUMO 
orbitals depends on the edge configuration, i.e., a- and 
f3 states can be non-degenerate or degenerate depend- 
ing on number of the carbon rings along the zigzag and 
armchair edges^^. The degeneracy reappears for large 
structures, such as n > 8, m > 7 (see notation in Fig. [T] 
(a)). The increase of the graphene size leads to disap- 
perance of the confinement effect and as a result closing 
of the gap. Thus, for n = 4 and m = 5 the gap is ~0.5 
eV and already for n — 6 and m = 7 the gap is sup- 
pressed to ~0.19 eV. The influence of the confinement 
effect on the graphene gap has been already confirmed 
experimentally^^. 

However, the state of D2h symmetry is not the ground 
state for nanoscale graphene. Graphene, optimized with 
C2v symmetry, where the mirror plane of symmetry is 
parallel to the armchair edges, has a total energy lower 
than that for the D2h symmetry. For the C2v symme- 
try, the HOMO and LUMO orbitals are characterized by 
the a— and /3-spin states localized on the opposite zigzag 
edges. Because the carbon atoms at the opposite zigzag 
edges belong to different sublattices, such spin distribu- 
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a) graphene b) doped graphene 



FIG. 1: (color online). The spin density distribution: (a) 
for nanoscale graphene optimized with the D2h point-group 
symmetry, and (b) for the case when one edge is doped by 
nitrogen, where the highest symmetry is the planar symme- 
try. Different colors indicate the a- (light) and /3-spin (dark) 
states. The spin density is plotted for isovalues of ±0.01 e/A^. 
The n and m are introduced to identify the the number of the 
carbon rings along the zigzag edge (n) and along the armchair 
edge (m). 



tion breaks sublattice symmetry and opens a gap (^ 1.63 
eV for n = 4 and ni = 5). The size of the gap is compa- 
rable to that found for nanoribbons^S. The large gap of 
nanoscale graphene obtained here is a result of significant 
contribution of the confinement effect, as the nanoscale 
graphene is confined in all directions. 

The localization of the a— and /3-spin states be- 
longing to HOMO and LUMO at the opposite zigzag 
edges is important for the application of graphene in 
spintronic o^'^'^^i'^° . However, the C2v symmetry state is 
a highly metastable state. Its total energy is comparable 
to that of D2h symmetry with a difference of ^ —0.5 eV 
for small structures such as n = 4 and m — 5, but the dif- 
ference decreases exponentially down to ~ —0.02 eV with 
increasing the structure size up to n > 6 and to > 7, that 
has good agreement with earlier workii, and disappears 
when n > 8 and m > 8 . The competition of the C2v 
state with CI symmetry, which is not constrained to have 
spin ordering along the zigzag edge, is even more crucial 
because of almost identical magnitude of their total en- 
ergy. However, we found that the distortion or dissimi- 
larity induced along a single zigzag edge not only breaks 
the sublattice symmetry of the graphene, but can control 
the spin ordering, thereby stabilizing the ground state of 
the C2v symmetry. The highest possible symmetry of the 
doped graphene is lowered from D2h to C2v symmetry as 
a result of the edge dissimilarity. The spin density distri- 
bution for the nanoscale graphene with one zigzag edge 
doped by nitrogen is presented in Fig. [1] (b). The local- 
ized states along the zigzag edges are formed by unpaired 
electrons belonging to the natural non-bonded orbitals, 
which participate in formation of HOMO and LUMO or- 
bitals. The a- and /3-spin states of HOMO and LUMO 
orbitals are spatially separated, i.e. localized at opposite 
zigzag edges. The (HOMO-1) and (LUMO-hl) orbitals 
usually correspond to the surface states, redistributed 
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over the entire graphene structure. The surface states 
are important for conductivity of graphene in a trans- 
verse electric field, because the charge transfer between 
the spatially separated HOMO and LUMO orbitals may 
occurs through participation of the surface states. The 
electron density distribution for the edge states and the 
surface states is presented in Fig. [51 The slight difference 
between the a- and /3-spin surface states (HOMO-1) is 
due to doping of the left zigzag edge. The a- and /3-spin 
states remain spatially separated with increasing struc- 
ture size. 



III. HALF-METALLICITY OF GRAPHENE 

The edge dissimilarity allows us to explore the required 
properties, such as the semiconductor- type bandgap and 
localization of a- and /3-spin states at opposite zigzag 
edges. Moreover, the spatial separation of the a- and 
the /3~spin states resulted from doping of the single zigzag 
edge is stable in comparison to the water adsorption^^. 
Doping of a single edge shifts the band energies of the 
orbitals which are strongly locahzed at this edge. Such 
a shift provides an opportunity to obtain another useful 
property which is important for spintronics - the half- 
metallicity of graphene. For the HOMO or LUMO or- 
bitals, which are shown to be localized at the edges, dop- 
ing can create a strong non-degeneracy of the a- and 
/3-spin states, because these states are spatially sepa- 
rated and localized at the opposite edges. Moreover, the 
HOMOq and LUMO^ orbitals are localized at one edge, 
while HOMO^ and LUMOa at the other. If doping in- 
creases the bandgap A„ for the a-spin state, then the 
bandgap for the /3-spin state, in contrast will be re- 
duced, and vice versa. Therefore, doping induces the spin 
gap asymmetry in graphene. Materials exhibiting asym- 
metric gaps for the a- and /3-spin states where one gap 
is of semiconductor type while the other is an insulator, 
are known as half-semiconductor materials, but if one of 
them is metallic, the system is half-metallic. Therefore, 
by choosing the right doping we can achieve a stable half- 
metallicity in graphene which will be an important step 
forward for applications in spintronics. 

We have investigated the transformation of the elec- 
tronic structure of nanoscale graphene due to the induced 
edge dissimilarities. The results are schematically pre- 
sented in Fig. [21 For nanoscale graphene with the D2h 
symmetry, a small bandgap occurs due to the confine- 
ment effect. The HOMO and LUMO orbitals in this case 
are localized at the zigzag edges, but their electron den- 
sity is equally redistributed over both edges (see Fig. [31 
(a)). Termination of the left zigzag edge by hydrogen 
(see Fig. [31 (b)) opens a gap as a result of breaking of the 
sublattice symmetry, thereby lowering D2h symmetry to 
a stable ground state of C2v symmetry. The hydrogena- 
tion leads to saturation of the dangling a bonds at the 
terminated edge but does not significantly change the 
energy of the HOMOq, and LUMO/j states localized at 



this edge. The resulting non-degeneracy of the a- and 
/3-spin states is not large, and the HOMO-LUMO gap 
of the a-spin state (A^ =1.8 eV) is almost identical to 
that of the /3-spin state (A^=2.1 eV). The doping of the 
left zigzag edge by nitrogen (see Fig. [31 (c)) shifts down 
the orbital energies of the HOMOq. and LUMO^ states 
localized at the doped edge and results in a strong non- 
degeneracy of the orbitals. This leads to a slight enhance- 
ment of the HOMO-LUMO gap for the a-spin state up 
to Aq,=2.2 eV, but significantly decreases the HOMO- 
LUMO gap for the /3-spin state down to A^ =0.8 eV. 
The length of the nitrogen-carbon bond at the edges is 
found to be c?Ar_c'=l-35 A, which is similar to the carbon- 
carbon bonds dc_c'=1.39 A. Similar results are obtained 
for phosphorus impurities, where the gaps are Aa=2.0 eV 
and A^=0.9 eV. Phosphorus is, however, less useful be- 
cause of the large phosphorus-carbon bond (dp„c=l-78 
A) which can lead to destruction of the lattice. We have 
also investigated the possibility to dope the single zigzag 
edge of the nanoscale graphene by other impurities, such 
as oxygen and boron, but they are not as effective as 
nitrogen. The oxygen doping leads to strong delocaliza- 
tion of the electron density of the orbitals localized at 
the edges. The doping by boron leads to even more trou- 
bles due to the long boron-carbon bonds at the edges 
((iB_c=l-42 A) and shifting of the states localized at 
the edges from the HOMO-LUMO gap deeper into the 
conduction and valence bands. 

For the nanoscale graphene structure investigated in 
this work, the spin asymmetry is achieved but bandgap 
magnitude for a- and /3-spin states corresponds to the 
half-semiconductor behavior (Aq,=2.2 eV,A_a =0.8 eV). 
Increasing the size of the graphene results in a decrease of 
both the Aq and A^ gaps due to the diminishing of the 
confinement effect. Therefore, for graphene structures 
doped by nitrogen or phosphorus of size n > 6 and m > 7, 
the A^ gap is closer to metallic type. Thus, for n — 6 
and m = 7 the gap for a-spin state is suppressed down 
to 1.13 eV while for /3-spin state down to 0.19 eV, which 
corresponds to the half-metallic behavior of graphene. An 
external electric field applied between the zigzag edges 
has been shown^i2i^i2£ to shift the band of graphene 
with spatially separated and degenerated a- and /3-spin 
states. The electric field shifts the bands in such a way 
that for the a-spin the HOMO and LUMO levels move 
closer to each other in the energy scale, while for /3-spin 
they move apart. At a certain electric field Aq van- 
ishes, thereby creating a metallic behavior of graphene. 
If +Ec electric field can close the bandgap for the a-spin 
state, the —Ec leads to bandgap disappearance for the /3- 
spin state. Therefore, the current voltage characteristic 
of such a structure will be symmetrical because the A^ 
equals A^ and for both spin states the switch from the 
semiconductor behavior to metallic occurs at the same 
critical electric field ±Ec- The advantage of graphene 
with spin gap asymmetry, i.e. different A^ and Ap gaps, 
found in this work is the different values of the criti- 
cal electric field required to close these gaps, such that 
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FIG. 2: (color online). Spin polarizations in nanoscale graphene where the left edge is doped by nitrogen. Different colors 
correspond to different signs of the molecular orbital lobes. The electron densities are plotted for isovalues of ±0.02 e/A^: (a) a- 
state of HOMO (^^homo = -6.04 eV) (b) /3-state of HOMO (Shomo = -5.43 eV) (c) a-state of (H0M0-1)(£homo-i = -6.38 
eV) (d) /3-state of (HOMO-1) (-Ehomo-i = -6.43 eV). The HOMO and LUMO are found to be localized at the single zigzag 
edges (edge states), while (HOMO-1) and (LUMO+1) - delocalized over the entire graphene structure (surface states). Bottom 
pictures show the representation of the localized and surface states. 
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FIG. 3: (color online). Schematic diagrams showing the distribution of the edge states and surface states in the energy scale 
and over the graphene structure (see the bottom pictures in Fig. [5] for pictorial description of the states) . The structures at 
the bottom demonstrate the spin distribution with isovalues of ±0.01 e/A^. (a) nanoscale graphene optimized with the D2h 
symmetry, (b) with left edge terminated by hydrogen and (c) with the left zigzag edge doped by nitrogen. For the localized 
states the energy levels (HOMO and LUMO) show density distribution (schematic), particularly delocalization of the orbitals 
between the two edges if the D2h symmetry is preserved, and their localization on the zigzag edges when sublattice symmetry 
is broken and C2v symmetry becomes to be highest possible symmetry. The surface states ((HOMO-1) and (LUMO±l)) are 
delocalized over the entire graphene structure (see for example Fig. [2] (c,d)). 



1 i?c(/3) l<l -E'c(a) I when < Aq,. Therefore, this struc- 
ture will be characterized by the spin-polarized current 
and by a non-symmetric current-voltage characteristics 
as for a semiconductor diode, when the current flow in 
one direction is preferable to the other. 



IV. DOPING OF GRAPHENE 



We have also investigated the influence of impurities 
on the electronic structure of graphene in the case when 
they are not embedded at the zigzag edges. Replacing 
carbon atoms by nitrogen atoms in a graphene lattice re- 
sults in the appearance of impurity levels inside of both 
the Aq, and A^ gaps. The energy diagram of localiza- 
tion of the molecular orbitals for the doped graphene is 
presented in Fig. 3] As we mentioned earlier, in pure 
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FIG. 4: (color online). Schematic diagram showing the dis- 
tribution of the edge states (LUMO, HOMO-1) and states 
localized by the dopant in the middle of the graphene struc- 
ture (HOMO) (see the bottom pictures in Fig. [5] for pictorial 
description of the states). The HOMOa and HOMO/3 are 
extra impurity levels that appear due to the doping and re- 
places the occupied orbital localized on the left carbon edge 
by shifting it dipper into the valence band. The inset picture 
(in brackets) demonstrates the electron density distribution 
for the HOMOc, with isovalues of ±0.01 e/A^. 



graphene the HOMO and LUMO orbitals are localized 
at the zigzag edges. The applied doping creates one ex- 
tra occupied orbital (HOMO) which is localized at the 
embedded nitrogen atoms and located above the occu- 
pied orbital belonging to edges, which becomes HOMO- 
1. The NBO analysis has shown that this extra orbital 
is formed by the unpaired sp electron localized on each 
nitrogen atoms. This reduces the HOMO-LUMO gap 
(Aq,=1.1 eV and Ap=Q.7 eV) while preserving the spin 
asymmetry. 



In an applied in-plane electric field the charge trans- 
fer occurs between the orbitals localized on the opposite 
zigzag edges, i.e., in our case such transfer occurs between 
the HOMO-1 and LUMO orbitals, which is a multi-step 
process with participation of HOMO. Because the gap is 
decreased and each nitrogen atom adds an unpaired elec- 
tron into the system due to the doping, the conductivity 
of graphene would be significantly enhanced, and can be 
controlled as it is done in semiconductor devices. 



V. CONCLUSION 

We have investigated the possibility to control the elec- 
tronic and magnetic properties of nanoscale graphene. 
We found that if pure graphene can be characterized by 
a small bandgap and no spin ordering at the zigzag edges 
the dissimilarity of the edges induced by doping impuri- 
ties lowers the highest possible symmetry to C2v, which is 
characterized by the spin ordering along the zigzag edges 
and their antiparallel alignment between opposite zigzag 
edges. Moreover, impurities embedded at a single zigzag 
edge shifts in the energy scale the molecular orbitals lo- 
calized at this edge, thereby decreasing the bandgap for 
one spin channel and increasing the other. Under these 
conditions, the half- metallic behavior can be achieved. 
Nitrogen doping in the middle of the graphene surface 
is found to have the prospect for application in nano- 
electronics due to the appearance of the occupied impu- 
rity levels in the bandgap. The impurity level results 
in a decrease of the bandgap of ^ 2.0 eV by one half 
and contains unpaired electrons, which should lead to an 
enhancement of the conductivity. Therefore, both the 
conductivity of the nanoscale graphene and its magnetic 
properties can be controlled by the impurities. 
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